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SIMULATION  AND  PERFORMANCE  OF  BRUSHLESS  DC  MOTOR  ACTUATORS 


SUMMARY 

The  simulation  model  for  a  Brushless  D.C.  Motor  and  the  associated 
commutation  power  conditioner  transistor  model  are  presented.  The  necessary 
conditions  for  maximum  power  output  while  operating  at  steady-state  speed  and 
sinusoidally  distributed  air-gap  flux  are  developed. 

Comparisons  of  simulated  model  with  the  measured  performance  of  a  typical 
motor  are  done  both  on  time  response  waveforms  and  on  average  performance 
characteristics.  These  preliminary  results  indicate  good  agreement.  Plans 
for  model  improvement  and  testing  of  a  motor-driven  positioning  device  for 
model  evaluation  are  outlined. 
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INITIAL  DISTRIBUTION  LIST 


-  LIST  OF  SYMBOLS 


Bm  Viscous  friction  coefficient  of  the  motor  (oz-in/rad  per  sec) 

Bi  Viscous  friction  coefficient  of  the  load  (oz-in/rad  per  sec) 

EgA  Back  EMF  Phase  A  (volts) 

EBB  Back  EMF  of  phase  B  (volts) 

Egc  Back  EMF  of  phase  C  (volts) 

IA  Phase  A  current  (amperes) 

Ig  Phase  B  current  (amperes) 

Ic  Phase  C  current  (amperes) 

Im  Power  supply  current  (amperes) 

Jm  Moment  of  inertia  of  motor  shaft  (oz-in-sec2) 

Moment  of  Inertia  of  load  shaft  (oz-in-sec2) 
kb  Back  EMF  constant  (volts/rad  per  sec) 
kt  Torque  constant  (oz-in/ampere) 

L  Inductance  of  the  stator  winding  (Henrys) 

Iab»  jCA*  IBC  Loop  currents  (amperes) 
kadj  Air-gap  flux  adjustment  factor 

kwi  Magnitude  of  speed  versus  torque  curve  slope  (rad  per  sec/oz-in) 
Ra  Resistance  of  the  stator  windings  (ohms) 

Rs  Power  supply  interval  resistance  (ohms) 

$A ,  ,  <t>c  Per  phase  air-gap  flux  (webers) 

Rq  'On'  resistance  of  the  power  transistor  (Ohms) 

Tm  Motor  restraining  torque  (oz-in) 

Tj  Load  restraining  torque  (oz-in) 

TAt  Tg,  Tq  Per  phase  developed  torque  (oz-in) 
wm  Speed  of  the  motor  shaft  (rad/sec) 


OR  Gear  ratio 


8  Angular  displacement  (radians) 

VAf  Vg,  Vq  Phase  terminal  voltage  (volts) 

Vg  Power  supply  terminal  voltage  (volts) 

RPSA,  RPSB,  RPSC  Position  sensors  logic  signals 
PWM  Pulse  width  modulation 


INTRODUCTION 


Recent  improvements  in  rare-earth  magnetic  materials  for  use  in  brushless 
dc  motors  have  allowed  reconsideration  of  electro-mechanical  actuator  systems 
for  applications  requiring  very  high  ratios  of  torque-to-inertia .  The 
investigation  discussed  herein  has  been  concerned  with  characterizing  mathe¬ 
matically  the  dynamical  features  of  a  missile  fin  actuation  system,  from  the 
input  to  the  brushless  dc  motor  to  the  output  shaft  of  the  mechanical  actuator. 
The  physical  model  is  based  upon  an  existing  prototype  actuator  currently 
under  evaluation  at  the  Naval  Weapons  Center,  China  Lake,  California. 

In  general,  brushless  dc  motors  produce  torque  through  the  interaction  of 
a  magnetic  field  generated  by  a  permanent  magnet  rotor  and  a  dc  generated 
magnetic  field  in  the  stator.  The  rotating  permanent  ma?net  eliminates  the 
rotating  armature  and  the  mechanical  wear  normally  associated  with  brushes. 
These  motors  fall  in  the  class  of  Permanent  Magnet  Motors  and  enjoy  certain 
advantages  over  wound-field  types  such  as: 

"...Linear  torque-speed  characterise cs ,  high  stall  (accelerating) 
torque,  no  need  for  electric  power  to  generate  the  magnetic  flux  and 
a  smaller  frame  and  lighter  motor  for  a  given  output  power"  [l]. 

Additionally,  the  brushless  dc  motor  is  characterized  by: 

"...controllability  over  a  wide  range  of  speeds,  capable  of  rapid 
acceleration  and  deceleration,  convenient  control  of  shaft  speed  and 
position,  no  mechanical  wear  problem  due  to  commutation  and  better  heat 
dissipation  arrangement"  [)]. 

The  fundamental  requirement  of  an  electro-mechanical  actuator  -control 
system  is  to  provide  torque  to  an  output  shaft,  sense  the  position  of  the 
shaft  and  adjust  the  torque  to  balance  the  load  when  the  desired  position  is 
reached.  This  must  be  accomplished  with  a  minimum  of  frictional  resistance 
and  delays  associated  with  the  inertia  of  the  mechanical  components.  Effects 
of  viscous,  static  and  coulomb  friction,  together  with  the  torque  required  to 
accelerate  the  mechanical  components  of  the  system,  lead  to  a  reduction  in 
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torque  available  at  the  output  shaft  and  an  associated  reduction  in  system 
performance. 

One  approach  to  the  analysis  of  the  electro-mechanical  actuator  system 
has  been  to  divide  the  system  into  two  sequential  problem  areas.  The  first 
deals  with  the  dynamic  analysis  of  the  brushless  dc  motor  and  development  of 
the  transfer  function  necessary  to  duplicate  actual  steady-state  and  transient 
performance.  The  second  area  deals  with  modeling  the  mechanical  system 
elements,  taking  as  input  the  dc  motor  shaft  angular  acceleration  predicted  by 
the  motor  analysis.  The  mechanical  system  must  be  modeled  considering  the 
effects  of  friction  and  inertia  and  translating  the  rotational  motion  of  the 
brushless  dc  motor  shaft  to  the  output  shaft  of  the  actuator  for  application 
to  missile  maneuvering  control.  The  results  obtained  from  a  study  that  has 
placed  primary  emphasis  upon  the  latter  problem  area  -  the  modeling  of  the 
mechanical  drive-train  leading  to  the  fin  shaft  was  presented  in  Ref.  3  and  4. 

This  report  deals  with  the  first  problem  area  and  presents  an  overview  of 
the  work  done  by  MacMillan  [4].  In  Ref,  4,  MacMillan  developed  a  model  for 
the  output  circuitry  of  a  transistorized  power  conditioner  that  provided  the 
required  motor  commutation.  The  model  developed  in  Ref.  4  used  an  ideal  (zero 
Impedance)  power  supply.  This  report  extends  the  model  to  include  internal 
power  supply  resistance  and  also  identifies  areas  for  further  improvement  in 
the  system  modeling. 

•  * 

The  next  section  of  this  report  presents  a  brief  description  of  the 
system  followed  by  a  development  of  the  model  and  an  analysis  of  the  simulated 
motor  performance.  Results  of  the  simulated  system  are  then  compared  to  the 
measured  response  from  a  typical  commercially  available  motor  and  recommen¬ 
dations  for  further  improvements  in  the  model  are  outlined. 


-SYSTEM  DESCRIPTION 

GENERAL 

The  system  is  viewed  as  a  position  control  device  to  maintain  an  output 
angle  under  an  applied  hinge  moment  due  to  aerodynamic  forces  on  a  fin  or 
aileron.  The  motor  is  a  permanent  magnet  dc  motor  with  feedback  in  the  form 
of  back  emf  proportional  to  the  angular  velocity  of  the  motor.  The  block 
diagram  of  the  dc  motor  and  power  conditioner  is  shown  in  Fig.  1  [4],  The 
mechanical  actuator  and  drive  train,  as  currently  envisioned,  introduce 
various  inertial  and  damping  loads  together  with  an  aerodynamic  force  and  its 
associated  fin  hinge  moment  that  must  be  overcome  to  produce  output  motion. 
An  operational  block  diagram  of  the  load  torque  is  shown  in  Fig.  2  where  the 
hinge  moment  and  motor  shaft  angular  acceleration  are  viewed  as  inputs  to  the 
drive  train  [2].  Figure  3  is  a  schematic  of  the  drive  train  which,  as 
presently  constituted,  includes  the  motor  shaft  (leadscrew),  ball  screw 
assembly,  and  the  crank  which  is  keyed  to  the  output  (fin)  3haft.  Inertial 
loads  are  considered  individually  within  three  major  subdivisions  of  the 
actuator;  the  output  shaft  to  crank,  crank  to  ball  screw,  and  ball  screw  to 
leadscrew. 

A  detailed  analysis  and  development  of  the  mechanical  model  with  the 
associated  assumptions  are  contained  in  Reference  2.  As  stated  in  the 
introduction,  this  report  reviews  the  development  of  the  motor  model  reported 
in  Reference  4.  For  convenience  in  the  development  of  the  model,  MacMillan 
considered  the  load  seen  by  the  motor  to  be  caused  by  the  inertia  and  kinetic 
(viscous)  friction  of  the  motor  shaft  as  described  by  Thomas  in  Reference  5. 

The  schematic  diagram  of  the  power  supply,  power  conditioner  and  motor¬ 
load  is  shown  in  Figure  4  where  it  is  noted  that  the  power  conditioner  logic 
and  driver  circuits  are  modeled  as  ideal  on-off  devices  with  zero  time  delay. 
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Figure  I  Motor  and  Power  Conditioner 

Simulation  Diagraa 


Additionally,  MacMillan  considered  the  power  supply  to  be  ideal  (zero  interval 
resistance)  to  further  reduce  the  complexity  and  simplify  the  development  of 
the  power  conditioner  model.  However,  in  this  report,  the  model  used  for  the 
power  supply  includes  the  internal  resistance,  Ra.  The  effect  of  including  R3 
in  the  requirements  and  performance  of  the  system  model  is  presented  in  the 
section  on  Analysis  of  Simulated  Motor  Performance.  For  additional  informa¬ 
tion  on  the  development  of  the  D.C.  motor  model  refer  to  Reference  4. 

DEVELOPMENT  OF  SYSTEM  MODEL 

The  basic  simulation  by  Thomas  used  a  single  power  supply  and  super¬ 
imposed  the  phase  currents  to  produce  the  motor  torque.  Motor  drive  was  then 
realized  by  multiplying  averaged  armature  current  by  a  torque  factor. 

In  the  model  developed  by  MacMillan,  the  supply  is  considered  to  be  a 
split  supply  of  equal  voltages.  Armature  current  is  assigned  a  positive  sign 
if  it  flows  in  the  positive  direction  (i.e.  into  the  motor).  The  use  of  a 
split  power  supply  allowed  MacMillan  to  apply  circuit  reduction  techniques. 
Because  of  the  symmetry  that  resulted  from  the  split  power  supply  approach, 
the  complex  3-phase  bridge-type  circuit  simplified  into  a  two-window  network. 
The  development  of  the  reduced  circuit  begins  with  the  definitions  of  loop 
currents  as  shown  in  Figure  5  where  the  power  supply  voltage  of  2V  is 
represented  as  V+  and  V-.  The  simulated  model  equivalent  circuit  is  shown  in 
Fig.  6  with  node  N  defined  as  the  mid-point  of  the  power  supply  and  no'de  0  to 
be  the  mid-point  of  the  motor  windings  and  neither  of  these  nodes  are 
considered  to  be  at  ground  potential  (the  ground  is  identified  as  the 
Reference  node).  Note  that  the  power  supply,  VIN  -  VIF  -  VIB  -  2V  and  that 
all  the  network  variables  are  defined  using  the  CSMP  model  variable 
definitions  [6].  Given  the  assumption  that  the  network  is  balanced,  that  is, 
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all  transistors,  diodes  and  field  windings  are  approximated  as  being  alike  in 
characteristics,  it  is  possible  to  apply  the  circuit  reduction  technique  of 
Thevenin  and  obtain  the  simple  2-window  network  shown  in  Figure  7  [4],  In 
Figure  7,  the  Thevenin  equivalent  voltages  and  currents  are  defined  in  terms 
of  the  CSMP  model  variables. 

The  back  emf  voltages  of  each  phase  ( VEMFA ,  VEMFB  and  VEMFC )  are 
developed  and  presented  in  Reference  4  in  which  it  is  shown  that  these 
voltages  are  summed  2  at  a  time  in  proper  sequence  to  compute  the  loop 
currents  and  resulting  phase  torque.  The  total  developed  torque  is  then  the 
sum  of  the  2-active  winding  phase  torques  over  the  proper  60  degree  of 
mechanical  angle.  Figure  8  shows  the  back  emf  voltages  generated  across  2 
windings  taken  2  at  a  time.  The  60  degree  of  mechanical  angle  that  produces 
torque  is  shown  by  the  logic  unit  output  levels  of  the  position  transducer 
RPSA ,  RPSB  and  RPSC  in  Figure  8. 

Reference  4  presents  a  detailed  development  of  the  power  conditioner 
model  that  includes  the  assumptions  for  the  switching  transistor  dynamics  used 
in  the  power  conditioner  model  as  well  as  the  development  of  the  harmonic 
air-gap  flux  used  in  the  motor  model. 

ANALYSIS  OF  SIMULATED  MOTOR  PERFORMANCE 

The  proper  evaluation  of  a  system  model  requires,  quite  naturally,  a 
comparison  of  the  actual  system  response  for  the  same  type  of  input  function 
used  for  the  model.  For  the  brushless  D.C.  motor,  it  is  assumed  that 
measurements  can  be  obtained  for  the  following  voltages  and  currents: 

-  Power  supply  terminal  voltage;  V3g 

-  Voltage  across  two  of  the  three  windings;  V^,  VbC,  and  Vca 
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-  Voltage  from  one  winding  terminal  to  ground;  Va»  Vb  and  Vc 

-  Power  supply  current;  Im 

-  Phase  current;  Ima,  Imb,  Imc 

Before  taking  data  from  the  simulated  model,  it  is  necessary  to  set  all 
constants  and  parameters  to  the  nominal  values  given  by  the  supplier  of  the 
motor.  Care  must  be  taken  to  properly  interpret  the  values  given  by  the 
manufacturer. 

In  particular,  the  value  of  the  Back  EMF  coefficient,  kb,  needs  special 
attention.  If  the  motor  data  sheet  indicates  that  this  coefficient  is 
obtained  by  measurement  of  peak  to  peak  voltage  across  two  windings,  then  in 
the  model  that  uses  distributed  sinusoidal  air-gap  flux,  an  adjustment  factor, 
kadj  is  required  to  insure  that  the  Back  EMF  values  are  in  agreement.  It 
follows  that  the  torque  coefficient,  kt,  must  also  be  adjusted  by  this  3ame 
scale  factor  [7].  The  value  of  kadj  can  be  computed  from  the  measured  no-load 
data  for  the  motor  current,  speed  and  applied  voltage  by  using  the  average 
characteristics  balance  equation, 

^m  “  (ysg  "  kbwnl)/(Ra  *  Rq) 
where  kb  -  kbm  kadJ 

kbm  »  manufacturer  supplied  Back  EMF  constant. 

Thus  the  adjustment  factor  becomes 

kadj  *  ^vsg  ~  *m(Ra  +  Rq)]/(kbmwnl) 

•  • 

For  example,  given  a  no-load  speed  of  3060  RPM  with  a  terminal  voltage  of 
30  volts  and  a  no-load  current  of  0.30  amperes,  kadj  -  0.825  (given  a  total 
series  resistance  of  1.47  ohms).  The  value  of  kadj  raust  of  course  be  less 
than  unity,  otherwise  the  motor  current  would  be  negative  in  value  which 
implies  that  generator  rather  than  motor  action  is  taking  place. 


Before  an  attempt  is  .made  to  close  the  loop  on  the  motor-load  to  form  a 
positioning  device,  it  is  important  to  operate  the  system  as  a  velocity  device 
and  obtain  measurement  of  the  motor  average  performance  characteristics. 
Typical  motor  average  performance  curves  are  speed  versus  torque,  motor 
current  versus  torque  and  output  power  versus  torque  from  no-load  torque  to 
near  stall  torque  conditions.  Figure  9  shows  these  three  curves  for  the  model 
with  torque  given  in  oz-in  units.  Note  in  particular  the  straight  line 
characteristics  for  motor  speed  and  current  and  that  the  power  output  is  of 
quadratic  form.  The  condition  of  constant  input  voltage  is  used  when 
gathering  data  for  these  curves,  and  it  follows  that  if  speed  is  a  straight- 
line  versus  torque,  then  current  also  has  a  straight-line  when  plotted  against 
torque.  This  can  be  shown  a3  follows: 

wm  *  wnl  ~  kwlT 
where  wm  »  motor  speed 

wni  -  no-load  speed 

kwl  -  magnitude  of  slope  of  the  speed  vs  torque  curve 

T  -  load  torque 

The  current,  using  average  values,  can  be  written  as 

^m  ”  (vsg  ~  kbwm)/(Ra  +  Rq) 
where  v3g  -  power  supply  terminal  voltage 

kb  ■  Back  EMF  coefficient 

•• 

Ra  »  winding  resistance 

Rq  -  transistor  "ON"  resistance 
then  by  substitution, 

i-m  ”  ^vsg  ~  kb  wnl)/(Ra  +  Rq^  +  t(kbkw^)/(Ra  +  Rq)]T 
where  the  first  term  on  the  right  hand  side  of  the  equation  represents  the 


no-load  current. 


OUTPUT  POWER  (WATTS)  VS  LOAD  TORQUE  (OZ-IN) 


The  power  output  versus  torque  curve  in  Figure  9  reaches  a  peak  value 
between  no-load  torque  and  stall  torque.  The  analytical  development  for  this 
quadratic  curve  is  presented  in  Appendix  A  where  it  is  shown  that  the 
condition  for  peak  power  output,  P0,  is  that  the  torque  must  be  of  value 

T  ■  (2kbwnl  ~  vsg ^ ^ ( 2kbkwl ^ 
to  produce  maximum  power  output 

PQ  *  ktvsg^(  **Bskb^ 

Since  the  ratio  k^/k^  is  a  constant,  increasing  peak  power  output  depends 
upon  decreasing  the  power  supply  interval  resistance,  R3,  or  increasing  the 
power  supply  voltage,  vag  as  one  would  logically  expect.  It  must  be 
understood  that  the  assumption  used  in  the  above  analysis  wa3  that  the  air- 
gap  flux  was  a  constant  average  value.  For  the  actual  motor  the  air-gap  flux 
is  distributed  sinusoidally  and  current  in  the  windings  produces  a  field  that 
may  cause  distortion  in  the  air-gap  flux.  At  large  current  values  that  occur 
at  and  above  peak  power  output,  the  field  distortion  may  result  in  an  increase 
or  a  reduction  of  output  power  for  a  given  load  torque. 

Another  motor  coefficient  to  consider  carefully  is  the  no-load  viscous 
friction  constant,  Bm.  The  value  of  Bm  is  often  not  included  in  the  motor 
data  sheet  and  the  probable  reason  is  that  its  value  will  depend  somewhat  on 
the  manner  in  which  the  motor  is  attached  to  the  system  load.  Bm  can  be 
calculated  from  given  no-load  data  as  follows: 

Bm  ”  (kt  kadj  W/wnl 
where  im  is  the  no-load  current. 

RESULTS  AND  CONCLUSIONS 


The  preliminary  results  indicate  that  the  balanced  bridge  circuit 
approach  used  in  the  development  of  the  Brushless  DC  motor  model  produces  good 


agreement  with  measured  motor  characteristics  as  indicated  below.  Further 
evaluations  of  the  model  as  well  as  improvements  and  additions  to  the  model 
will  be  conducted  in  the  near  future  using  data  currently  being  gathered  from 
a  prototype  of  a  fin  positioning  actuator. 

A  comparison  of  model  and  motor  produced  Back  EMF  voltage  is  shown  in 
Figure  10.  The  upper  curve  is  the  simulated  motor  waveform  and  the  lower 
curve  is  the  measured  Back  EMF.  The  model  adjustment  factor,  katjj  was  set  to 
0.63  value  by  trial  and  error  until  peak  to  peak  voltages  were  in  agreement. 
Another  verification  of  the  model  is  obtained  by  comparison  of  Figure  11  and 
Figure  8.  Both  figures  show  waveforms  of  Back  EMF  across  2  windings  taken  2 
at  a  time  and  also  show  the  timing  waveforms  for  the  Position  Sensor  Devices 
(RPSA,  RPSB  and  RPSC)  for  counter  clockwise  rotation.  The  waveforms  agree  in 
both  phasing  and  in  form. 

The  steady-state  performance  curves  of  a  typical  motor  are  given  in 
Figure  12  where  the  motor  load  torque  is  given  in  lb  in  units.  These  curves 
were  produced  with  a  constant  terminal  voltage  of  30  volts  and  agree  in  form 
with  the  curves  of  the  model  as  given  in  Figure  9.  Peak  power  output  occurs 
at  load  torque  of  approximately  95  oz-in  for  both  motor  and  model. 

Additional  validation  of  the  model  is  shown  in  Figure  13  where  the 
current  in  Phase  C  for  the  model  (upper  curve)  is  in  close  agreement  both  in 
form  and  in  phase  with  the  same  current  for  the  motor  (lower  curve). 
Additional  model  results  are  shown  in  Appendix  B  where  typical  input  voltage 
step  response  waveforms  are  presented. 

FURTHER  RESEARCH 

Improvements  in  the  model  are  required  with  regard  to  the  triggering  on 
of  the  diodes  used  to  protect  the  transistors  from  excessive  reverse  currents 


PHASE  C  CURRENT  -  AMPS 


Figure  13  Phace  C  Current:  Top  Curve-Simulated  Mddel, 
Lower  Curve-Typical  Motor 
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daring  the  switch  off  time  for  the  transistors.  MacMillan  in  Ref.  was  able 
to  provide  this  action  for  the  condition  that  the  power  supply  interval 
resistance,  Rs  is  zero  value.  When  Rs  is  not  zero,  additional  programming  is 
required  to  perform  the  required  trigger  action. 

The  Pulse  Width  Modulation  (PWM)  used  by  Askinas  in  Ref.  8  must  be  added 
to  the  present  model,  however  Askinas  used  a  separate  additional  transistor  to 
perform  PWM.  The  Power  Conditioner  to  be  used  in  the  Test  Stand  for 
measurement  of  actuator  performance  accomplishes  PWM  by  using  the  lower  3et  of 
the  commutating  power  transistors.  Therefore  some  modification  of  Askinas* 
CSMP  coding  may  be  required  to  accurately  perform  the  PWM  speed  control. 

Extensions  of  the  model  will  Include  a  position  sensing  device  and  a 
tachometer  as  well  as  the  nonlinear  loading  as  developed  by  Franklin  in  Ref  9. 
Test  data  results  will  be  compared  to  model  in  both  transient  (step)  response 
and  in  frequency  (Bode)  response  and  the  results  will  be  reported  in  a  Summary 
Report  for  FY1986. 
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the  rate  of  change  and  set  it  equal  to  zero.  Thus, 


k.k  ,v 
t  wl  st 


dT L  rA  +  Rq 


2k. k.  k  ,w  , 
t  b  wl  nl 

rA  +  Rq 


2k.  k.k  .2 

t  b  wl  _  0 

RA  ♦  Rq  *L 


Solving  for  load  torque,  TL  yields 

*  “nl^^wl  “  Vsg^^^g^wl^ 
tL  *  (2kfcWni  -  vSg)/(2k{jkwi) 

and  substitution  into  the  power  output  equation  produces  the  conditions  for 
peak  power  output  as  follows: 

P0  "  k^Vgg/(Ra*Rq )  [ wn^-kw^ (2k{jWn^-vag )/(2k{jkw| )  ] 

~  k^k^/ (Ra+Rq)  [wnl“kwl(2ktjWni-Vgg)/2kt)kwi )  ]2 
After  expansion  and  cancellation  of  like  terms,  the  peak  power  output  reduces 


P0  ■  kfcVgg/[4k(j(Ra+Rq)  ] 

As  a  numerical  example  consider  the  following  values  of  a  typical 
Brushless  DC  Motor  operating  at  a  steady-state  speed  of  3060  RPM  (320.4 
rad/sec)  with  a  constant  terminal  voltage  of  30  volts. 
kt  -  15.9  oz-in/ampere 
kb  «  0.112  volts/rad  per  sec 
kadj  *  0-825 

Ra  ■  1.37  ohms,  Rq  -  0.10  ohms 
kwi  «  1.67  rad  per  sec/oz-in 
Thus  for  this  example  motor, 

PQ  -  05. 9)(0.825)(30)2(7.062x10"3)/[4(. 112X0.825X1  .47) 

PQ  *  153.5  watts 
and  occurs  at  load  torque 

To  -  [(320.4)/(1 .67)  -  (30.)/(2)(.’12)( .825) ( 1.67)] 

T^  *  94.65  oz-in  ■  5.92  lb-in 


APPENDIX  A 


Power  Output  of  a  Brushless  DC  Motor 

The  mechanical  power  developed  by  a  Brushless  DC  Motor  is  given  by  the 
equation, 

P  -  (7.062x10~3)  Twjj  watts 

where  T-Tl+Tm,  the  sum  of  the  load  torque,  and  the  motor  restraining 
torque,  T^  in  oz-in  units  and  wm  is  the  motor  speed  in  rad/sec.  units.  At  no 
load  TL  -  0  and  at  stall  torque  wm  -  0,  thus  power  reaches  a  peak  value  in  the 
region  between  no  load  and  stall  torque  conditions. 

The  evaluation  of  steady-state  peak  power  can  be  obtained  analytically  in 
terms  of  the  load  torque  by  substitution  as  shown  below. 

Given  that  T-Tl  that  is  Tl>>T^  and  T-ktim  with  wM-wNL  where  kt  is  the 
motor  torque  coefficient  and  kw^  is  the  slope  of  the  speed-torque  curve. 

Then  the  power  output  is 
po  ■  ^t^m  "  kwlTL.) 

and  since  steady-state  speed  conditions  are  assumed,  the  current  is  given  by 

v  -  k.w 
,  sg  b  m 
1"’ 


where  v3g  is  the  voltage  at  the  terminals  of  the  power  supply,  k^  is  the  Back 
EMF  coefficient,  Ra  is  the  armature  resistance  of  2  phases  in  series  and  Rq  is 
the  sum  of  2  power  transistors  forward  'conduction  resistance.  By 
substitution, 


kt  v  K.  K. 

P  .  t  -3-S  fM  -  |<  T  )  -  — — -  (w 

°  Ra  +  Rn  1  nl  Kwl 1 L '  Ra  +  Rn  ^  nl 


kt> 
A  +  nq 


k  ,T.  ) 
wl  L 


and  to  find  the  conditions  for  peak  power  output,  it  is  necessary  to  compute 


A1 


The  curve  of  power  output  versus  load  torque  for  this  example  motor  is 
shown  in  Figure  9.  It  must  be  noted  that  the  above  analysis  is  valid  for  the 
assumption  that  motor  current  produced  fields  do  not  have  a  measurable  effect 
on  the  magnitude  and  phase  of  the  flux  generated  by  the  rotating  permanent 
magnets  and  also  that  a  constant  average  flux  exists  in  the  air-gap  between 


rotor  and  stator. 


_•  “  Jl 


APPENDIX  B 

CSMP  Listing 

The  system  described  by  the  schematic  diagram  in  Figure  4  is  programmed 
in  the  CSMP  language  as  shown  below.  The  computer  used  was  an  IBM  370  located 
in  the  W.  R.  Church  Computer  Center  at  the  Naval  Postgraduate  School. 
Following  the  program  listing  is  a  User's  Guide  for  the  recommended  procedure 
in  setting  the  constants  and  parameters  for  validation  with  measured  motor 
performance.  Also  included  are  typical  output  waveforms  for  a  step  input  of 


supply  voltage. 


FILE:  brtlEA 


csmp 


A  1 


//JBBTFO0 1  JCd  (C 16  3,C J2  2)  .  '  KACH1I 1  A-  12  •  ,cl  ALS  =  G  ,  MSGLEV£L=  ( 1  ,  1) 

//•MAIN  CSG=N?JV«1.C1fcSi<ilNiI=SC5 

//•FORMAT  PR,DCNAflE  =  ,D££i=I  C  q  A  L 

//•FORMAT  EB,DLNAHE  =  5YSVECl£i,LESI  =  ICCAL 

//CSMPSTEP  EXEC  CSHPVD V ,PCF 5l R £  =  M  \  C  I F 5 

//X.SYSIh  CC  * 


PATCHES  FRCH  IOC  INS1ILLEC 

VERSION  ELEVEN  —  SIXTH  REVISION  I  NCCFPCFA1ING  1HESV5  LIU. 

THIS  VERSION  1NCCBPCBAIES  REVERSING  CC  t  BUI  AT  ION  AI  FEV1IH. 

CURRENTS  ARE  NCI  Sl’FZFFCSEE.  ELUX  •  CUREEM  IS  CCHPCIEC  FCR  EACH  PEASE 
AN G  THE  RESULTING  TCFCUES  AFE  SUMMED. 

THE  PURPOSE  OF  VERSION  FIVE  IS  1C  TREAT  THE  FLUX  AS  VARYING 

ACCORDING  TO  THE  SUB  CE  A  FUNDAMENTAL  SINUSCIC 

AND  ITS  FIFTH  HARMONIC  AS  EXPLAINED  IN  CHAPTER  FIVE. jTECBAS) 

THE  TOTAL  FLUX  IS  APPFCXIHAIEC  AS  THE  AIGEERAIC  SUN  CE  THE  FLUX 
DEVELOPED  IN  1-0  WINDINGS  AI  A  T IB  E  IC  CALCULATE  THE  TCSCUE 
COEFFICIENT  AND  LACK  EMF  CC EFFICIENT. 

THIS  MCDEL  SIMULATES  MAGNETIC  FLUX  AS  A  SINUSOID  AND  A  F IF  It 
HARMONIC  AND  SIMULATES  JIOLE  CCMHUIATICN  AS  NELL  AS  THE  SNITCHING 
LCGIC  AND  TRAN SIS  ICR  DYNAMICS  CE  A  fcUUSHLESS  DC  MCTCR  PCNER 
CONDITIONER.  THE  NINCINGS  ARE  NCI  TFEAIEL  I N C EP E NE E Nil Y 
AND  THEIR  TRANSIENT  ELECTRICAL  INTERRELATION  IS  SIMULATED, 

THE  CCNTFIfiUIICNS  1C  DEVELOPED  1CFCUE  IRE 
TBEATED  AS  SUP E P PC S A E I E . 

THE  PROGRAM  HAS  EEEN  MODIFIED  1C  INCLUDE  PCNER  SUPPLY  RESISTANCE ,  ES. 
NITH  THE  ADDITION  CF  PS,  THE  1FIGGE  PING  VCITAGE  FCE  THE  TRANSISTOR 
PROTECTING  DIODES  HAS  StEN  CHANGED  TC  INSURE  PROPER  CPEFATICN. 


INI1IA1 

* 

CONSTANT  KTK  =  15. Hb,  l !.  =  0.01465.  EL  =  0.00,  JL  =  O.C,  N  *  1.0,  ... 
JM  =  0.CU1,  Kell  =  J.II^J,  PI  =  2.l41592c5,  K  K  2  *  3.5SC,KADJ  =  .6j,  ... 
THADV=  3.5,  ?S  =  €.22,  ICIMC  »  C.C1 

PARAMETER  LA  =  .COOd,  FA  =  0.6£5C,  111  =  C.C 
PARAMETER  VSAT  =  .4,  RSAT*  .C5,  cCCI  =  1.CE+4 
PARAME1ER  IT  IMS  =  1.0E-6,  EEV1IH  =  C.13 

PARAMETER  V D 1 D  =  U. , V D2 D  =  C . , VD 3D  =  0. . V C4C=0 . , V D5 C  =  C. . VD6C  =  C  . 

PARAMATER  VSGFEF  =  JC.C,  KINT  =  1CCCC.3,  K  li  A  K  P  =  18250. C 


KTM  --  MFG  SUPPLIED  BC10R  IORCUE  CCNSIANT  (CZ-IN/AMP) 

KT  —  PHASE  ICPCUE  CONSTANT  IC2-IK/AMF) 

KBH  —  MFG  SUPPLIED  MC10F.  HACK  EMF  CuNSTANI  ( VC  L 1/ R  A  D/ S  EC  ) 

KE  —  PHASE  SACK  EMF  CCNSIANT  (VCIT/FAD/S) 

KADJ  --  ADJUSTMENT  FACTOR  FCR  MFG  SUPFIIEl  SCTCF  CONSTANTS 
BA  —  FHAiE  RESISTANCE  OF  THE  MC1CF  (CHM) 

LA  —  PHASE  INDUCiANCE  CF  THE  MC1CF  (HENRIES) 

dH  —  VISCOUS  FFICTICN  COEFFICIENT  CF  THE  MCICH  (C  l - 1 N/ F A  C/S ) 

BL  —  VISCOUS  cSIoTICN  COEFFICIENT  CF  THE  LCAL  (C Z- 1 N/ E AC/ S) 

BLP  —  VISCOUS  FRICTION  COEFFICIENT  CF  LOAD  IhilU  F.EEUC1ICN  GEARS 

e  —  Total  viscous  fpicticn  cf  ire  motor  system  (see  helcw#nc-scfi) 

JM  --  INERTIA  CF  THE  SCICr  (CZ-1N/S-S) 

JL  —  INERTIA  CE  THE  LOAD  . 

JLP  —  INERTIA  OF  THE  LOAD  IHFU  REDUCTION  GEARS 

J  --  TCIAL  INERTIA  CF  THF  MCTCF  SlSTEX  {SEE  EELCW,  NC-SCFT  S5CIICN) 

A2  =  J/E  —  THE  MECHANICAL  TIME  CCNSIANT  CE  THE  MCTCR  (SEC) 

ABTAU  =  LA/  ( ?  A  ♦  F  EC  A  D )  —  IHE  ELECIFICAI  TIME  CONSTANT  CF  THE  MOTOR 
AND  LnlVc  TrANSISTCRS  FCF  CUFFEM  PATH  A-E 
BCTAU  =  LA/ <F A+FECtC)  —  THE  Ei.fcCloI.cAL  TIME  CONSTANT  CF  THE  MCTCR 
AND  DRIVE  iRANSiSIO^S  FCR  oUFRFM  PATH  ?-C 
CATAJ  =  LA/  ( 5  A  ♦  R  E  w  C  - )  --  Hit  ELECT.-IOAl  TIME  CCNSIANT  OF  IhE  MCTCR 

AND  D-.UVE  IF  ANSI  SI  OPS  FCP  CUFi-rM  PATH  C-A 
V I F  *  \/t  SUPPLY  VOL  I  AG  r. ,  VIN,  MEASUltJ  r?oM  ♦  To  MIC  £CSIIICN,N 
VIB  =  1/2  SliFrLY  VOLTAGE,  VIN,  MEASURED  FHCM  -  IC  MID  ECSIT1CN,N 
BC»  --  EQUIVALENT  olRCJII  CilSiANCt  CE  >.  (») 

CfSX?  -  ■■’O  .«£."<  ilAL  MODEL  PC  i  T  R  A » S  1ST  TC «  EETWEFN  3ATUSAIICN  AND  CUICFF 


FILE:  D.i'ihA 


*  RSAT  = 

*  RCU1  3 

*  TIIflE 

*  IH6SI 

* 

*  THADV 

*  THCCN 

*  8S - 


tl.-cJII  FSSISIANCI  Cf  1BANSI3ICR  SA1U5A1ICN 
CiFC  JII  'VrGISlANCE  CE  IBANSISICa  CUICEF 


ECUiVALEGi  C  I .« C  IJ 1 1  FE3I21ANC1 
E-CIV.ii.ENI  CiFC  JII  RESISTANCE 
--  StaliChlNI  II  !*  E  CF  TSANSISICK 


_  1  j.  v-  n  x  »  j  .irt  v-r 

—  THIS  IS  IHDEG  n  HICn  IS  SE1  1C  22FC  AT  STA5I.  SESEI 
AI  3d 0  DEG  AND  CAE  ALSC  E£  fJi£C  AS  SI.NUIATICN  SICE 

—  THETA  ADVANCE  FC  3  CCHSOIATTCN  (CCi.  IS  FCSIIIVE) 

—  COdMOIATION  CCNTHCL  ANGLE=1HESI»1KACV  IN  DEGREES 
INTEBNAL  EESISI ANCE  C i  TEE  SUPPIK  VC11AG5 


1C  2EBC 
ANGLE. 


NOSCE1 

Kl=KTd*K AC J 
KB=KBM*KACJ 
BLP  =  EL/ 

JLP  3  JL/  (N**.<) 

J  3  JM  ♦  JLP 

fc  =  ctl  ♦  ELF 

AI  =  LA  /  BA 

k2  -  J  /  t 

A3  3  LA/ (BA  ♦RSAI) 

B  S 1-  B  3/  2 . 0 

BS2=BS/2. 0 

VINIC  =  15.0/KINT 


ME1HCD  SlliE 
DYNAMIC 

VIF=  VlNbAF 
Via  3- VI fcd  AF 
VIN  =  V I F  -  VI 

*  VSGDEL  3  (VSGF 

*  VSGE8R  3  DEADS 

*  V INI  NT  -  IN'IGB 

*  VINHFG  =  BIN! 
VINHAF  3  VINEF 


VRAMP  1 
VRA ME2 
VINBF  1 
VINHF2 


(VSJfEF  -  VSG ) / 
DEADS? (  -3.5  ,  ♦ 
IN'IGBL  (VINIC,  V 
BIN!  ♦  VI  Mil 
VINEF1  -  VINFF2 
RAdPiC.O) 

BAN?  (J.Og  11 
NRAMr^VRAMr  1 
KRA.M  *VF  A.1F2 


-  VSG)/2.J 

-3.5,  *3.5,  VSGDEL) 

I  NIC .  VrCtnB) 


*  INTRODUCE  FINITE  iiANSixICN  1 1 .1 E  FCR  IfANSlSlOT-  SWITCH 

*  BE  INCCFPC RATING  HXf  CNEMIAL  RISE  AND  CECAE  IN1C  Si.  1  —  S  fc6 

*  giEXP  PROVIDES  A  REALISTIC  EXPONENTIAL  TFANSI1ION  EETBEEN  cUTCFE 

*  AND  SATCRATICi  -iilCH  UCiUDHS  SAIUBAIICN  DELAY  iihEN  fLI  IHFCGGfl 

*  THE  LIHIIEB 

w 1  EXE  =  NcALPL  (.75.  ilidE,  3.75-SM) 

22£XE  =  n  £  ALtL ( J .  7o ,  ill, ME,  J.75-S*$) 

Q3EXE  =  PEALPL  ( J. 76  ,  T  LIME,  3.75-s*J) 

C  4  E  X  £  -  itEAi.PI  (0.7 5  ,  1 1 1 M  E  ,  J.75-SN4) 

Q5EXP  3  uEALPL  ll. 76,  TUNE,  3.75-S*5 
goEXE  =  •iEAi.Pl  (G. 76,  HIM,  0.75-SX6) 

EEL23R  ;1EXP  =  0.1,  *2cX?=3.1,  Q3FX?3C.l,  «'MXI  =  G.l,  g5£XP=C.1,  s .fcEXi 
ABSEBR  1  £  X  P  =  0  .  1 ,  2^EXE  =  0.1,  v3EXt  =  C.1,  g4EXE=C.1,  g55XE  =  C.l,  *'cEX 
RC1  =  Lilli  (.-SAT,  RCC1,  2  .  *FC’J  I  ♦;  1  E  X  t) 

Rw^  =  L1.11I  (foal ,  NCCi,  2.  ♦FCJI+22E  (?) 

RQ3  =  LIMIT  ( B  S  A I  ,  BCCI,  2 . ♦ SC U  I  * C 2 E X i ) 
aC4  *  LI. Ill  (BSAT,  ECU!,  2  .  *  FC  J 1  *  C  4  £  X  ?) 

8g5  3  LIMIT  (BSAT  ,  RCU1,  *BCU1  *w  5F  Xt:J 
Rg6  3  LIMIT  (BSAT,  FCU1,  2  .  * rC U T ♦ C c 2 X ?) 

!  CICDE  MODEL 

.r  n  t  **  » i  h  r  t  t  *  /*  r 


LINEAR  CICDE  MODEL 
DIODE  LUuN-CN  VCIIAGES 
VD  ID  3  fcCCO.O  -  VAIND  ♦ 
V  CD  D  3  D  u  w.  C .  ♦  V  C  I  N  C  ~ 

VD3D  3  tCCC.C  -  VwInC  ♦ 
VC2D  3  fcCuO.J  ♦  VALID  - 
VD5D  3  t  3  C 0 . 3  -  VC  INC  ♦ 
VD4D  3  cCSJ.C  ♦  VLINC  - 
P  D  1  3  ECUS-  (Vine,  iv  S  A  L  , 
RD2  3  SC3SB  (VL^C,  cSmT, 
ECJ  3  ICNo  A  (  VD  2C  .  l-SAl. 


VCINC 

vein: 

V  A  I  Jt  C 
Vv.  INC 
v  i : :.  c 
v  Ait.: 

2  CUT, 
SCUT, 


FILE:  (MIS  A 


Ci  lr 


R Do  =  iCN5i(  (VbtCf 


:iCU,  FC  L  I ) 


1  HE  IHEV2NIN  r.  J  U  I V  A  L  E  N  ‘I  Cl  TC  (J 1 1  KEIJGD  IS 
VOLTAGE  !'ID-PClVi,N,  AS  RErEFENCE  AND  IHE 
APPLIED  IC  COMPUTE  IHE  PHASE  CIFFEMS. 


USEC  5UIH  IEt  SUPPLY 
LCCP  MEIriCC  CF  ANALYSIS 


THEVENIN  ECUIVAIENT 
VANTfi  =  VN  l  *  BECAS2/ 
VBNTH  =  V(i1*E£CBS2/ 
VC NTH  =  VNl*PEwCSi/ 


VCITAGES 


EECASI  +  IECASi)  ♦  V  N2  *  B  E  C  AS  1/  S  EG  AS  1  ♦  F  E  2  AS2 
BECHS1  ♦EECES2)  ♦  VN2*FtCES1/  1  EF<,ES1*EECES2 
SE<2CS1*P£CCSi)  ♦  VN2*FECCS1/  SEqCS  1  +  B  E  wCSZ; 


TEANSISICp/DILCE  E.UIVriLENI  6ES1SIANCE 
REQA1  =  RC  1  ♦  PC  1/  (Pi  HtiC  1) 

R E Q A 2  -  R  C  ^  *  RE2/(rC2*BD2) 


REQB1  =  RC3  +  R£ J/ + 
BEQB2  =  6^4  *  FLU/  ic^H  +  r.ZU 
BEQC1  =  *  RE5/(E*5*RD5 

B  EqC  2  =  R  1 6  ♦  r<C6/<ri_6  +  RDft 


r*E6/  (f^6*FD6] 


CON  POT  A1 1C  N  CF  IilEVENIN  FCU1VALENI  RESISTANCE 

REQAS1  =  EECA1  ♦  RSI 
SE^ASE  =  n  E  C  A  2  ♦  FS2 
REGBSl  =  BE«E1  ♦  BS 1 
REQBS2  =  RECc2  ♦  FS2 
HECCS1  =  fctjd  ♦  F S 1 
R£yCS2  =  BEQC2  ♦  RS 2 
THEVENIN  EQUIVALENT  BESIS1ANCE 
REQA  =  B  E  2  AS  1  ♦  F  EC  AS  2/ (R  EC  A  S  1  ♦  F  LC  AS2) 

RE JB  =  BtCESI  *  FEwfc£2/  ( F  E  v  dS  1  ♦  F  EC  t  S  2) 

FEQC  =  BECCS1  *  Br;CS,£/(RECCSl*F£QCS2) 

TOTAL  LCCP  RESISTANCE 
REQAE  =  2 . *R  A  ♦  BEqA  ♦  C  EC  E 
REyBC  -  2  •  *  a  A  ♦  FECI  ♦  F.  ECC 
RECCA  =  2.  ♦Rd  ♦  I  ECC  ♦  I? E C A 

LOCF  COBBER  I-  MAXIMUM  VALUE 
IA3  =  (VANId  -  V  c  N  T  fc ) / £  E  C  A  c 

IAB  A  =  MVIN  ♦  VEMEis  -  V  *  M  E  A )  /  (BS  ♦  2.*3A  ♦  BEC^i 

t  —  *  t;  "  »  it  t’  nn  (t  r  v  r  /■ 


IBC  =  ( V i:N  1 H  -  V C  N1  i. )  /  ?  E C c C 

IBCA  =  -  ( V I N  ♦  VEXES  -  VEMEC)/tFi 
ICA  =  (VCNIB  -  VANTf )/5ECCA 
IC  AA  =  -(VIN  ♦  VE1FC  *  v!*FA)/|£S 

TIME  CC NS  I  AN  IS 
AEIAU  =  <..*  LA/ NEC  AT 
BCTAU  =  2.*  LA/RECcC 
CATAU  =  *.♦  L  A/  1 1  „  C  A 

LOOP  CUFFENTS 

1MAB  =  FEALPL  (1.C,  AtTAU,lAL) 

1H3C  =  SEAL2L  ( C .  C , ec I A U, IcC) 

IflCA  =  B thLPI  ]C. 0,CAI AO, IC A) 

NET  LEC  (PhASZ)  CURFENIS 

IMA  =  IKAt  -  INCA 

IM3  =  IMSc  -  IMAB 

IMC  =  1«CA  -  IMEC 

IMCD  =  EERIV  (J. O.IMC) 

INCUR  =  TCIMC^IflCC 

IMCAC  =  hEALPL  (Q.O,TCIMC,ir.CL£) 

IMF  =  FEA1PL  (0.0.0. CC5.IM) 

VSGF  =  BEALPI  ( 0 . 5 , 0. C C ‘ , V  SC ) 

INDdCICR  l.iEUCEE  VCITAGES 
LI  ADI  =  DE«  U  {  2.  C,  I.'A) 

DIuDT  =  D  E R  I V  ( C  .  0 ,  I/.  L-) 

CICCI  =  LES;V(C.«J,  I«C) 

VAINE=L A*[ lAJi 
Vd1NC  =  LA*C  ltd 
VClNC  =  LA*£ic£  L 
PHASE  CBM'S  LA.  VILLAGE 


2.  *F  A 
2.  *R  A 


B  E  CC  1 
BECA1 


BEC  A  1) 
BECE2) 
FECC2) 


VCITAGES 
,  I  *  A  i 


g  a  v  -  n 

0,  I  - 


VILLAGE 


^  »  *  V  ^  A  - 


FILE:  EMU  ,1  CSMd  A  1 


VAOLV=IMA*  (ft  A) 

VEQL  V  =  IM  d  *  (BA) 

VCOLV=lKC*  (BA) 

*  PHASE  VCLIaGE 

V AO  -  ISA* (BA) ♦LA*CIAC1*VEEFA 
VBO  =  ISa* (flA) ♦LA*EIECI*VEftFc 
VCO  =  IMC* (HA) ♦IA*EICEI* VEMrC 

* 

*  CCA  PUT  AIICH  OF  POHEB  SUPPLY  TEFMliVAL  VCLIAGES  (GNC  FEF) 

*  WHEBE  VIN  IS  SUPPLY  VOLTAGE  ANC  IM  IS  oUPtLY  CUkFESI 

VSG  =  VIN  -  I«*FS 

VAG  =  VSG  *  F.ZC  A2/  (r  Es,  A1  ♦  HECA2) 

VBG  =  VSG  *  Eii*c2/ (FECcI  ♦  F  EC  1 2) 

VCG  =  VSG  *  FE\,C2/(FECC1  ♦  3  EC  C  2) 

VAC  =  VAG  -  VCG 

*  VOLTAGES  ACROSS  1BAN5ISICSS  ANC  EICEES 
VQD 1  =  VSG  -  VAG 

VCD2  =  VAG 
VQD3  =  VSG  -  V£G 
VCD4  =  VEG 
VQD5  =  VSG  -  VCG 
VCD6  ^  VCG 


P V AG  =  LIMIT  (-60.  ,  EC.,  VAG 
PVBG  =  L 121 II  -60.,  fcC . ,  VEG 
P  VCG  =  L 1 .1 1 1  -60.,  6  C  •  ,  VCG 

P  VAC  =  rVAG-FVCG 
F  VBA  =  EVEG-EVAG 
PVCB  =  PVCG-PV2G 


igi  = 

FCNSW ( 

[3  W  1 , 

.003,  . CC2, 

IMA) 

IC2  = 

r  C  N  S  » 

5  W  2 , 

.  o  0  3  ,  .  0  G  d  , 

IMA) 

IC3  = 

FCNSW  j 

Si3, 

. C  C  3  ,  .  302, 

IMF1 

IQ4  = 

i  CNS  ’«  { 

5  W  4  , 

.00 3,  .003, 

I  ME 

igs  = 

ECNSW  J 

S  W  5  , 

.002,  .003, 

I  MC 

IC6  = 
PC  1  = 
rg2  = 
PC  3  = 
SC4  = 
FQ5  = 

FCNSW  I 

V  L  C  1  * 

V  CO  2  * 

V  E  C  3  * 
VG-Z4  * 
VC* 5  * 

S  «b  . 

I  v* 

1:1 

IC4 

10- 

.003,  .003, 

IMC 

PC6  = 
PQTOT 

vegb  < 
*  ?£1 

INC 
♦  PC2 

♦  P 0 3  ♦  PC 

4  ♦  i 

PC  t 


*  THE  PHASE  ANGLES  LOCK  cCNEUSEC  BUT  IhE  PHASE  FEIATICNSE1F  IS  FIGHT 

BEilFA  -(3.+2iN(2.*IhciA*(11*?l/t))  ♦.5^*SIN(1C.*THllA*J1*El/6))) 
BEflE  b  =  (3.* SIN  (2.  *1H£I  A*  (7<?I/c)  )  ♦.S^  +  SIN  ( 10.  ♦THETA*  (5*El/o) ) ) 
EEMFC  = (3.*S1N (2. *ihE i A*  j*Pl/6)  ♦ . j9*SI S ( 1C. ♦THE! A*  »9*PI/fc) ) ) 

*  NORMALIZE  LEG  EME 

VEMFA  =  t  E  M  F  A  *  KB/KKd  *  w  M 

VE1FE  =  cEMFE  ♦  KE/KK3  *  W  M 

VEMFC  =  d  E  MFC  *  K  C  /  is  K  3  *  WM 

VEMFAC  =  VEMEA-VEMEc 
VEMFcA  =  VEMfc-VSMEA 
VEMFCE  =  VEilEC-VEMFt 


*  ICRCUE  IC  MAKE  MCTCF  TURN 
TA  =  IMA  *  K1  *  EEMEA/KKJ 
TB  =  I  ME  *  KT  *  EEMFE/KK j 
TC  =  IPC  *  K1  ♦  EEflf C/KK3 
TM  =  TA  ♦  TE  ♦  TC 

TflM  =  Ki  ♦  IM 
TBM  =  EM  *  » M 
VKHbt  =  Kc*«K 

*  1CBCCE  LNIEGPAL  VALUE 

*  1 M I N I  =  I  NIG  E  L  (C.  0  ,  IM) 

^  1L  =  TLL*3IFi  (1.0 £-i) 

T  N  2  =  TNI  *  (1. 3/t) 

WM  =  ft  E  A  L  £  i  (d. Ct Az,  I;.<) 
WMtNT  =  IMJrL  (  j.  :,»'•) 


B5 


:  Vt-julJi 


FILE:  EM  If  A 


Li.1t  A 1 


*  WH=320. 4 '4 i 5 

nCJRPM  =  WK  *  1.3  2  ./FI)  **. 

fcMRPMF  =  -MtifVrt 
THS1A  =  LNIGRLjO.O.hfll 
1HOEG  =  LHETA  *  (18C.C/PI) 

* 

*  COMMUTATION  ADVANCE  HESE 
1HCON  =  THRST+IHAOV 

* 

PWfifll  =  WM  *  18  *  . OC  7  3625 
PWHM  a  EM  *  . 0 C7 Oc  7  5  *  wM  **2 
PHE  =  EL  *  .  ) 0  7 'jc 7 6  *  »M**z 
IL  =  VIE/ (2*6601) 

*  THIS  PROCEDURE  FRCVIEES  A  SIMPLE  MECHANISM  fCf  REVERSING  THE  MOTOR'S 

*  DIRECTION  WHILE  KEEPING  THE  LCAl  AS  AN  CEECSING  1CRC0E. 

PBOCEDUHE  IN1  =  fWCEfcE(iiB,TM,II) 

IF  (NB.LT.C.O)  GC  1C  1 
TNI  =  1H  -  II 
GC  TC  5 

1  INI  =  IB  ♦  I L 
5  CONTINUE 
ENDPRCCFCUBE 

*  THIS  £  5CC  E  DUE  E  RESETS  THE  VAFIAELE  THRil  TC  C  AFIFa  EVER*  360  IEGPEES 

*  OF  MECHANICAL  RCTAIICN.  IBIS  IS  FUNDAMENTAL  TC  THE  SIfillAlICN  CF  ALL 

*  SWITCH  I NG  AND  POSITION  SENSING  AcTICN. 

PROCEDURE  1HRSI=RESET (1FCEG) 

THESI  =  AMCC  (T 6 D E J ,  2oJ.) 

IF  (1HRST  .LI.  0.0)  TURSI  =  TiiFSI  ♦  360. 

ENDEHCCEDURE 


*  THIS  PROCEDURE  WAS  IAKE3  r  PC  1  IhE  NEXI  PROCEDURE 

*  SENSCHS  SHvJULD  GENERATE  A  SINGLE  VALUFE  FUNCTION 

*  REGARDLESS  OF  DIRECTION  CF  aCIAIICN 
PROCEDURE  St1,sE^.SE3  =  aALL  (I I: CCS) 

IF  (I  HOC  N  .  (>  &  .  130.)  GC  TO  45 
GC  TC  46 

4b  THCCN  =  IHCCN  -  1  ci  L*  - 
46  CONTINUE 


IF 

[THl.CN 

.  GE. 

U  .  .  A  N  D  . 

THCCN 

.  IT. 

jC. 

GC 

IC 

10 

IF 

THCCN 

•  ii  £  • 

tGm  .  A  N  L  . 

THCCN  .IT. 

c  0 . 

GC 

IL 

1  1 

IF 

THCCN 

•  biL* 

"  C  C  N 

.  LI. 

50. 

G C 

:c 

12 

IF 

IhCCN 

•  Cl  C.  m 

90.. A  SC. 

IHCCN 

.11. 

1*0. 

GC 

*c 

13 

IF 

THCCN 

•  CJ  C  • 

IHCCN  .IT. 

15  0. 

GC 

1L 

14 

THCCN 

.GE. 

15j..ANC. 

l  ii  C  C  !i  .LI. 

Id  0. 

GC 

TC 

15 

5£  2  =  0 . 
SE3  =  1. 
GOTO  2 C 

11  S  E 1  =  1. 

S  E  2  -  C  . 
5  E2  =  0 . 


12 


13 


14 


1b 


20 


*C 

1. 

1. 

0. 

-cO 

C. 


GCTC 
SE  1  = 

SE2  = 

SE3  = 
GCTC 
Stl  • 

SE  2  =  1. 
SE  3  =  0. 
GOTC  20 
SE  1  =  0. 
S  E  2  -  1 . 
SE  3  =  1. 
G  C  I  C  £.  v 
SE  1  =  0. 
S  E  2  =  C  . 
SE  3  =  1. 
CC NT  I  VUE 


cICAUSE  THE  HAIL 
CF  JrCSIIICN 


FILE:  Bill  HA 


Lilli 


A  1 


ENDFRCCECC5E 

*  THIS  PF.OCEDU?  E  SIMULATES  GENERAL  C  CMMUT  AT IC  N  , 

*  DETERMINES  IHE  A  LGE  E  F  A IC  SUM  CF 

*  THE  VAFIAELE  FLUX  CCttECN  ENTS  F  C  5  USE  IN  CCfPUIING  THE  HCICF'S 

*  APPROXIMATE  BACK  EMF  (BEMFT)  .  II  ALSO  S'JEIBACTS  THE  GZNEBA1ED 

*  VOLTAGE  FROM  THE  PBCEER  SUPPL*  VOLTAGE. 

*  Sil  1HEU  SW6 

*  SIHULATE  POWER  TRANSAS IGF  IRIGJE3S  BEING  ENERGIZED  CR  SWITCHED  OFF 

*  THCON  IS  THE  V1RIAEIE  THROUGH  WHICH  1  HE  SWITCHING  LOGIC 

*  IS  IBPLEHENTED.  REVTIM  IS  IKE  I1.1F  AT  WHICH  CLOCKWISE  CCHMUIAIICN 

*  BEGINS. 


PROCEDURE  SWl,SW2,SVJ,SW4,5W5,SW6,E£FFI,VN1,Vt*2  =  CC(1fl<IieE,... 
REVTIH .IHCCN.BECFA, c EHF E , dEHFC. VEPf A, VFSfB , VEMFC) 

If  (TIME  .GT.REVTia  )  GO  1C  1 0 C 
THCON 1  =  AMCC  (THCON ,  IdG.) 

If  (THCON 1  .LI.  C.C)  THCON  1  =  THtCNl  ♦  130. 

If  (THCON 1  .GE.  U..AND.  THCON  1  .LT.  30.)  GO  IC  50 


U..AND.  THCON  1  .LT.  30.)  GO  IC  50 


IF  (THCCN 1  .GE.  30..ANC.  IHCCN1  .IT.  6C.  GO  IC  51 


IF  THCCN1  .GE.  bC..AND.  THCCN  1  .LI.  SC. 
If  THCCN1  .GE.  9C..ANC.  THCCN  1  .11.120. 
IF  THCCN  1  .GE.120..ANE.  THCCN1  .L1.15C. 
IF  THCCN  1  .GE.  15u..ANC.  THCCN  1  .LT.IdC. 
*  CLOCKWISE  COrtMOIAIlCN 
100  CCNTINOE 

IF  (THCCN 1  .GE.  0..A.VC.  IHCGN1  . LT .  30.] 
IFlTHCCNl  .GE.  30.. AND.  THCCN  1  .LT.  fcO. 


LI.  SC.  GC  1C  52 
II. 120.  GC  1C  53 
11.15C.  GC  TO  54 
LT.IdC.  GC  TC  55 


IF  (THCCN  1  .  GE . 
If  (THCCN  1  .GE. 


IF  (T  HCC  N 1  . GE.12C..AND 
IFlTHCCNl  .GE.  150. .AND 


tO. .AND.  THCON  1  .IT.  SC. 

90. . AND.  THCGN1  .LI. 120. 
2C..AND.  IHCCN1  .11.150. 

50. . AND.  THCCN 1  .LI.  130. 


GC  TC  53 
GC  TC  54 
GC  TC  55 
GC  TC  50 
GC  TC  51 
GC  1C  52 


8 EM FI  =  EEMFC  -  EEMfc 
VN  1  *  V I F  -  V E K F C 
VN  2  =  VIE  -  VEMFB 

m«ihc 

GC  TO  fcC 
51  SW1  =  1. 

5W2  =  C. 

SW  3  =  G. 

SW4  =  1. 

SW  5  =  C. 

In  ^  —  0  • 

BEMFT  =  *  c E M F  A  -  £  £ M F  3 


VN 1  =  VI F 
VN2  =  VIE 
IH=  I  £1 A 
GC  TC  cC 
52  Stal  =  1. 


V  Z  M  F  A 

V  E  fl  F  E 


30 

5 J  SW  1 
SW  2 
SWJ 
3V4 
SW5 
SHfc 


BEHFI  =  EESFA  -  EEKFC 
VN 1  =  V I F  -  VEBF A 
VN2  =  VIB  -  VEBFC 
IH  =  ir.A 

30  IC  c C 

SW 1  =  c. 


UEKFT  =  EEMFE 


L  c  A  E  C 


1 5V  W.V  AWW  1  *- 1  >. !  I."  1/  I  ■  M 


FILE:  E.1TE-A 


CSM? 


A  1 


i 


V. 


VN  1  =  V  IF  -  VE.MrE 
V  N  2  =  Vic  -  VEMI-C 
IM=IRE 
GC  1C  toO 
sa  SHI  =  c. 

SH2  =  1. 

SH  3  =  1 . 

SHU  =  0. 

SHE  =  C. 

SHfc  -  0 • 

a  E  B  F  I  =  BEilEE  -  EERFA 
VN  1  =  VI F  -  VERrB 
VN 2  ^  Vic  -  VEBFA 
IH=IBE 
GO  10  cO 
55  SHI  =  0. 

SU  2  =  1. 

SH3  =  0. 

SHU  =  3. 

SVI 5  =  1. 

SHE  =  0. 

BEBF1  =  fcEBFC  -  EEHEA 
VN  1  =  VIF  -  VEflFC 
VN2  =  VIS  -  VEBFA 
IHSIHC 
60  CONTINUE 
ENDP30CEDUBE 

TEBMIN AL 

TITLE  BASIC  DC  BOTCE  SYSTER 

TIHEfi  FINIIM  =  .050.  CUTDE1  =  .0000275,  EHCEl  *  .00C275,... 

CEIRIN=  1 .  3£-  1C 
FINISH  THESI  =  UOO.C 

PRINT  nfl,  HMhEM,  lHESl,lNflHINI(lHINl#naiNI(IAVelAV,bflAV, 
VAIND,  V  fc  I ;»  D  ,  VCIND, V  N  1.V  S2,  V  AC  .  V  EC,  VCC  ,  C3  .  V  C&U  ,  V  £  D  5  ,  V£C  6  , 
VSG,  VAG, VEG, VCG, VD  ID  .VJ2D.VD3C. VCUL, VC5L  ,VECC,  ... 
IAB,IAEA,IdC,IECA£ICA,ICAA,IB,IECEG,  1 HCCR  , IhCGN  1.  ... 

I  BA,  I  St),  IRC,  EjBfVIN,VINHAF,irCAC,VSGF,lflfi,TEM,VK£bE 
LABEL  EdASE  C  CUBBSNl  -  ARES 
OUTPUT  TIRE,  I MC  ,  S  fc  2 
PAGE  XYELCI 
END 

&ESE1  £  &  I N  T 

LABEL  ECU  ER  SUEELY  CURRENT  TC  R.CICB  -  ARES 
CUTFUT  I  IRE,  I R  , S  «  4 
PAGE  XYELCT 
ESC 

LABEL  SCIur  SEE  EC  -  r.' E  S 
OUTPUT  IlnE, 

PAGE  XYELCT 
END 

LABEL  ECUS5  SJEILi  TERMINAL  VCI1AGE  -  VCLIS 
OUTPUT  TIRE,  V  3G , S  *  t 
PAGE  X'.PLCi 
E.NC 

LABEL  TCLAL  DEVELOPED  10R£(jE  -  C2-IN 
OUTPUT  TIDE,  1H,SS1 
PAGE  XYELCT 
EAC 

LABEL  EE-rILTEHEC  PHASE  C  ciJSRiM  -  ARES 
CUIPUT  ii.BE,  IRC  AC,  SI*  3 
PAGE  XYELCI 
END 

LABEL  VOLTAGE  ACROSS  A-C  TEFKIKALS  -  VCLTS 
OUTPUT  TIME,  V AC , SW 3 
PAGE  XYFLCI 
END 

LABEL  LE-FILIERE3  PC*  OP  SUPPLY  CURRENT  -ARES 
CUTPUT  TIME,  1  R.  F  ,  S  k  5 
PAGE  XYELCI 
c  N  C 

LAEEL  cACK  EMr  PHASE  C  -  VCLIS 
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FILE:  B n  I F A  CbJI? 


OUTPUI  lIHi,  V£EFC,Sii1 
PAGE  XiFICl 
£  N  £ 

LABEL  III  CLVELCFEC  ICtyUE  -  CZ-IN 
OUTPUT  TIME,  T f! fl  ,  S  i3 
PAGE  XYFLCT 
£  N  C 

LABEL  FCMEB  OUTCUT  -  fi A T T 5 
OUTPUT  TIKE,  PBB,SW5 
PAGE  XYFLCT 
END 

LABEL  LP-FILTEEEC  SUPPLY  TEH B INAL 
OUTPUT  TIKE,  VSGF,Si*1 
PAGE  XYELC1 
END 
STCE 

END JOB 

/* 

//DISSFCF.BYDIHS  CD  * 

DR  A  fc  =  1-END, 

flODIFY=  1-END  ( 3)  <S I Z £ =8 . , 2 . )  , 
BODIFY=2-ENC  (J)  (C  V £F F LOT  , CCE  N  E  F 
BOOLE Y =2- END  <J)  4GV2FFLOT,CChNEF 


VOLTAGE 


VOLTS 


<SIZ£=8. ,2.)  , 

(CVEFFLCT,CCENEF=0.,3.,SIZE=8.,3 

<CV£FELOT,CCHNEF=0.,3.,SIZ£=8.,3 


Users  Guide 


The  following  procedures  are  recommended  to  be  used  in  running  the  CSMP 
language  program  given  in  Appendix  B: 

1.  Obtain  the  motor,  load  and  power  supply/conditioner  data  from  the 
manufacturer  and  enter  the  values  in  INITIAL  section  at  the  beginning 
of  the  program. 

2.  Set  parameter  KRAMP  to  the  desired  input  voltage  V  as  follows: 

V  3 

example  V  -  30  volts,  KRAMP  -  xIO  -  15000.  Under  these  conditions, 

a  fast  terminated-ramp  input  results.  This  is  equivalent  to  a  step  input 
of  V  volts. 

3.  The  constant  KADJ  must  be  set  experimentally  to  produce  the  measured  motor 
Back  EMF  voltage.  To  do  this,  remove  the  asterisk  in  the  statement 

*  WM  -  320.4475 

where  320.4475  is  the  magnitude  of  the  driver  motor  speed  in  rad/sec.  (in 
this  case,  speed  -  3060  RPM).  Run  the  program  with  different  values  of 
KADJ  until  the  Back  EMF  across  2  windings  agree  in  magnitude  with  the 
measured  value  obtained  from  the  motor. 

4.  Set  N  to  the  value  of  the  speed  reduction  from  motor  shaft  to  load 
shaft. 

5.  Set  RS  to  the  value  of  the  power  supply  interval  resistance  and  LA  and 
RA  to-  the  per  phase  winding  values. 

6.  Set  RSAT  and  VSAT  values  for  the  power  transistor  "ON"  values  and  RCUT  to 
the  cut-off  resistance  value. 


7.  The  value  of  BM  should  ttfc  calculated  from  the  no-load  data  as  a 

preliminary  value  and  adjusted  experimentally  until  no-load  operation  of 
the  model  agrees  with  measured  no-load  values  for  the  motor.  Compute  BM 
as  follows: 

BM  -  ( KT ) ( KADJ ) ( IMNL ) / ( WMNL ) 

For  example 

BM  -  ( 1 5-89) ( .63) (0.3)/ (320. 4) 

BM  -  0.00937  oz-in/rad  per  sec. 

Notes:  1.  The  asterisks  in  column  1  of  the  program  statement  makes  the 
statement  a  comment  and  disregarded  by  the  CSMP  Translator.  These  statements 
are  the  result  of  program  development  and  are  included  for  future  application 
in  the  development  of  the  complete  actuator  model. 

2.  A  set  of  typical  output  waveforms  are  included  here  to  indicate 
the  plot  output  capabilities  of  the  program  which  are  in  the  TERMINAL  section 
of  the  CSMP  listing. 

3.  The  load  applied  to  the  motor  should  be  adjusted  over  the  range  of 
no-load  to  near  peak  load  and  the  performance  curves  of  Figure  12  plotted. 
The  CSMP  simulation  output  can  then  be  checked  against  these  results  by 
entering  BL  as  follows: 

BL  -  TL/WM  oz-in/rad  per  sec 

where  TL  is  the  load  torque  (oz-in)  WM  is  the  motor  speed  (rad/sec). 
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